A Spectral Investigation of Criticality and Crossover Effects in Two and Three Dimensions: Short Timescales with Small Systems in Minute Random Matrices

Random matrix theory, particularly using matrices akin to the Wishart ensemble, has proven successful in elucidating the thermodynamic characteristics of critical behavior in spin systems across varying interaction ranges. This paper explores the applicability of such methods in investigating critical phenomena and the crossover to tricritical points within the Blume–Capel model. Through an analysis of eigenvalue mean, dispersion, and extrema statistics, we demonstrate the efficacy of these spectral techniques in characterizing critical points in both two and three dimensions. Crucially, we propose a significant modification to this spectral approach, which emerges as a versatile tool for studying critical phenomena. Unlike traditional methods that eschew diagonalization, our method excels in handling short timescales and small system sizes, widening the scope of inquiry into critical behavior.


I. INTRODUCTION
The phenomenology of critical phenomena, encompassing phase transitions in diverse contexts, stands as a cornerstone within the framework of Statistical Mechanics theory.Initially conceived within the realm of many-body physics, it has evolved into a concept with far-reaching applications spanning disciplines such as economics, network theory, sociophysics, game theory, and numerous others [1][2][3][4][5].
In spin systems, a particularly effective approach to delve into critical phenomena involves conducting short-time dynamics studies.This method entails preparing systems with carefully chosen initial conditions and then analyzing temporal averages.Through this analysis, the manifestation of critical behavior in the system is revealed via power-law dynamics [6][7][8][9][10][11].
Utilizing time-dependent simulations to investigate these techniques presupposes the capability to compress temporal evolution while contending with inherently large systems, a challenge frequently encountered in equilibrium simulations.Moreover, equilibrium simulations grapple with the issue of critical slowing down.Nowadays simulations of two-dimensional systems typically encompass linear dimensions denoted by L, conservatively ranging in the order of hundreds or, with more generous allocations, extending into the thousands.These simulations leverage parallelization techniques and GPU applications to effectively handle computational demands.This suggests that exploring non-conventional methods is intriguing.Upon delving into the literature, the concept of random matrices reveals fascinating aspects.Originating from the intricate description of the distribution of energy levels in heavy nuclei proposed by E. Wigner [25][26][27][28], its connection with the thermodynamics of Coulomb gases primarily stems from the works of Dyson [29][30][31].
Alternatively, utilizing individual time series to construct appropriate matrices, known as Wishart matrices [32,33], offers a particularly intriguing avenue for delving into the statistical mechanics of spin systems and their idiosyncrasies.This approach is underscored by recent contributions [34][35][36][37], as these matrices encapsulate the time-correlations of specific random variables described by stochastic processes, which can be modeled by Langevin equations or simulated via Monte Carlo (MC) Markov chains.It's noteworthy to mention that while the integrability aspects of lattice spin systems have been explored using random matrices, the emphasis has traditionally been on direct examination of the Hamiltonians of such systems [38], rather than on correlation matrices as in the previously mentioned works.
In contrast to other methods found in the literature, the approach advocated by one of the authors of this paper in [36] diverges significantly by eschewing the use of matrices whose dimensions correspond directly to the number of lattice sites.Such an approach, while theoretically sound, can be computationally prohibitive due to the sheer scale of the matrices involved.Instead, the proposed method operates effectively with dimensions on the order of a few hundred, derived from various short-term evolution samples of the magnetization system.This strategy, as evidenced in our prior research [36,37], proves to be both computationally efficient and theoretically robust.
In this study, we aim to extend the applicability of random matrix methodology by investigating models that exhibit tricritical points, with a particular emphasis on the transition from critical points to the tricritical one.Through careful analysis of the spectra derived from Wishart-like matrices, which capture the correlations among diverse magnetization time series, we showcase the effectiveness of our approach in elucidating the phase diagram of the Blume-Capel (BC) model in both two and three dimensions.Throughout this paper, we present our findings, highlighting how this methodology adeptly reveals the intricate characteristics of complex systems like the BC model.
Furthermore, we demonstrate that this method, despite its efficiency and suitability for short times, operates effectively for small systems, contrastively to the standard MC method.
In the next section, we introduce the Wishart-like method alongside a concise overview of the established properties of the BC model.It's crucial to emphasize that our primary focus lies not on the intricacies of the model itself.Rather, our goal was to select the simplest model featuring a tricritical point in both two and three dimensions to validate our study.
In the following section (Section III), we unveil our principal findings, demonstrating the efficacy of our method in accurately pinpointing the critical points of both the two-dimensional and threedimensional versions of the model.Additionally, we illustrate how our approach adeptly captures the crossover effects in a spectral manner.Furthermore, a novel aspect explored in this study, not previously investigated in our prior works, is the successful application of extreme value statistics in establishing the criticality of such systems.
Lastly, we demonstrate the effectiveness of our method even for very small systems.We illustrate that although there may be a computational cost associated with diagonalizing matrices, this cost is offset by the advantages gained from shorter computation times and smaller system sizes.The paper concludes with a summary of key findings and conclusions.

II. RANDOM MATRICES, CRITICAL, AND TRICRITICAL POINTS IN THE BLUME-CAPEL MODEL
Tricritical points play a vital role in contemporary research within condensed matter theory and statistical mechanics.The pioneering discovery of the first tricritical point in He-3 and He-4 mixtures by R. B. Griffiths in 1970 marked a significant milestone [39].Subsequently, Griffiths, along with Blume and Emery, introduced the BEG (Blume-Emery-Griffiths) model in 1971, which provided a framework to replicate the thermodynamic behavior observed in these mixtures [40].This model, based on a spin-1 Ising model, has since become a cornerstone in the study of tricritical phenomena.
In its broader scope, the Hamiltonian can be expressed as: The initial term delineates the crystalline interaction among the spins, whereas the subsequent term denotes the multispin interaction among them.Finally, the third, fourth, and last terms represent the Ising interaction between spin pairs, anisotropic interaction, and the interaction of spins with an external magnetic field, respectively.
Indeed, an even more straightforward variation of this model, commonly referred to as the Blume-Capel (BC) model [41,42], is delineated by the Hamiltonian: where each spins can hold the values s i ∈ {−1, 0, +1}, it would be enough to present the existence of a tricritical point in both two and three dimensions, separating a critical line of the first order line.
The first term models the local interaction between the spins, with J > 0 representing the interaction strength and ⟨i, j⟩ indicating that the interaction occurs between nearest neighbor pairs of sites i and j.The parameter D is called the anisotropy field and is responsible for zero-field splitting, resulting in an increase in energy for s i = ±1 states even in the absence of an external magnetic field.Finally, the third term models the interaction of the system with an external magnetic field of intensity H, which we will assume is not present.
In this scenario, the model delineates a critical line (CL) culminating in a tricritical point (TCP).Subsequently, it exhibits a first-order line (FOL), as illustrated in Figure 1 for both the two-dimensional and three-dimensional versions of the BC model in the absence of an external magnetic field.
This figure is utilized only for pedagogical reasons in this work since all points used in this current work uses the points estimated in this line that were very didactically obtained by Butera and Pernici [43].Sure, a lot of estimates shown in such reference were re-obtained by the authors since there was a evolution of estimates with very different methods along of many years that include the a lot of good works (see for example [44][45][46][47]).Thus the question is if we can detect the critical line of the BC model for different values of D using Wishart-random matrices spectra and how the method responds to the crossover between CL and FOL intermediated by the TCP point.
Therefore, leveraging a framework established in prior research, this study demonstrates our capability to localize critical points within two and three-dimensional BC models, while also inves- tigating the existing crossover phenomena.Subsequently, the following section provides a concise overview of the random matrix methodology employed in our approach, along with the quantities slated for estimation within this method.
A. Wishart-like matrices and spin systems In our analysis, we introduce the magnetization matrix element m tj , which denotes the magnetization of the jth time series at the tth Monte Carlo (MC) step in a system comprising N = L d spins.For simplicity, we set d = 2, as it represents the minimal dimension for the manifestation of phase transitions in short-range interaction systems.In this context, t ranges from 1 to N M C , and j ranges from 1 to N sample , thereby constructing the magnetization matrix M with dimensions To delve into spectral properties, an intriguing approach is to shift our focus away from M and instead examine the square matrix of size N sample × N sample : where each element G ij of G is defined as t=1 m ti m tj , referred to as the Wishart matrix [32].To simplify computations, it is advantageous to transform the components of the matrix M using the transformed matrix M * , whose elements are expressed in terms of standard variables as follows: where This transformation facilitates subsequent analysis and calculations. Thereby: where Here it is crucial to expound upon a pivotal calculation that elucidates the application of these matrices in greater detail.We consider two distinct time evolution samples of the magnetization per spin denoted as m ti and m tj , where t = 1, ..., N M C .In this context: where σ t,j,k denotes the value of the k-th spin in the j-th evolution or run at time t.
We can establish the correlation between these two time series using the following definition: When T > T C , ⟨m i ⟩ t ≈ 0. This leads to: , and we can express the correlation coefficient (our matrix element of G) as: where σ i ≡ (σ i,1 , ..., σ i,N ) and σ j ≡ (σ j,1 , ..., σ j,N ).Thus g ij for T > T C is determined by: This metric assesses the relationship between the temporal averages of spatial correlations within both inter and intra-time series.By analyzing both spatial and temporal dimensions, it provides a compelling approach to delve into spin systems.
Thinking in the general case, when the variables m * ij are uncorrelated random variables, momentarily forgetting the context of these variables represent magnetization of spin systems, the eigenvalue density ρ(λ) of the matrix G * = 1 N M C M * T M * conforms to the well-known Marchenko-Pastur (MP) distribution [48].For our specific case, we express this distribution as: where . Undoubtedly, we expect that for T ≫ T c , the density of eigenvalues ρ exp (λ) obtained from computational simulations approaches ρ(λ) in Equation 3, but our method may not necessarily fit such a distribution perfectly due to residual autocorrelation.The interesting question is what happens when T ≈ T C .Moreover, we will utilize the density ρ exp (λ), obtained from computer simulations, to determine the critical parameters of spin models.
The moments of ρ exp (λ) are calculated as: where N bin is the number of bins of the histogram of ρ exp (λ).Thus, for T ≫ T c , we also expect λ k to approach: . However, beyond these limits, the behavior of λ k can provide thermodynamic information about spin models, as suggested by our previous works [36,37].In those works, we observed that monitoring ⟨λ⟩ and (∆λ) 2 = (λ − ⟨λ⟩) 2 as a function of T T C indicates a minimum of ⟨λ⟩ and an inflection point for (λ − ⟨λ⟩) 2 (or divergence of its derivative) occurs at T = T C .
Here, we will demonstrate that this method works effectively for the BC model in both two and three dimensions, particularly in identifying critical points and examining its response to the crossover phenomena between CL and FOL.

III. RESULTS
Now, we will present our main results.In the first subsection, we showcase the outcomes of our spectral method concerning the critical points of the BC model in both two and three dimensions.An essential aspect for the accurate numerical application of the method involves utilizing the histogram to compute the eigenvalue moments through numerical experiments, as per Equation 4, and directly calculating the numerical moments.It is crucial to emphasize this point for readers intending to apply the method, as we have confirmed that computing the averages directly does not yield the expected results presented here.

A. Critical points
We begin our results by displaying the histogram of eigenvalues.We choose D = 1 for both the 2D and 3D BC models to illustrate the density of eigenvalues obtained through the diagonalization of matrices G * .Figure 2 presents histograms for various temperatures.An evolution of the gap between the two eigenvalue bulks can be observed.Similar behavior is noted for the Ising model on two-dimensional lattices under mean-field approximations [36,37].
We can discern analogous behavior in the three-dimensional BC model, as illustrated in Figure 3.
It's crucial to highlight the distinctive trend of the eigenvalue gap narrowing around the critical temperature, along with the correspondence to the MP law for T > T C .While a perfect agreement is not expected as T approaches infinity due to the correlation matrix's construction, incorporating total magnetization and time series with inherent autocorrelation, it's important to note that this doesn't diminish the method's validity in any manner.However, it is necessary to utilize this density of states to effectively determine the critical parameter.This can be achieved by computing the moments of the density of states, specifically ⟨λ⟩ and (∆λ) 2 .In this regard, we observe the results for three different values of D. For the two-dimensional BC model (refer to Figure 4), we tested three values: D = 0, D = 1, and D = 1.75, employing the corresponding T C values estimated in [43] as a basis.Similarly, for the 3D BC model, we utilized D = 0, D = 1, and D = 2.2, as depicted in Figure 5.
We can observe a pronounced minimum in ⟨λ⟩ at T = T C in both the two-dimensional and three-dimensional versions of the BC model, which is related to the closing gap observed in Figures 2 and 3. Additionally, an inflection point seems to be observed for the variance exactly at T = T C in both versions of the model (in two and three dimensions), demonstrating that both spectral measures -the average and variance -are effective in exploring criticality.The inset plot displays the first derivative of the variance: indicating that the critical temperature is associated with a pronounced minimum (a negative value of significant magnitude), where t = T T C To better understand such behavior, we examine the second derivative: in the gap between two bulk eigenvalues is observed compared to the two-dimensional BC model (see Figure 2) and its plot is depicted in Figure 6 for both scenarios: the two-dimensional and three-dimensional BC models.
We notice that the critical temperature precisely aligns with the inflection point of the eigenvalue variance due to the condition ζ < 0 for T < T C and ζ > 0 for T > T C .Understanding the nature of this inflection point warrants further investigation, prompting a thorough discussion.In our work, we provide an in-depth analysis of this aspect in the appendix.
Thus, in this first subsection, we observed that critical points of the BC model are well captured by this spectral methodology in both versions of the model: two and three dimensions.We used different parameters based on fluctuations of the eigenvalues and their convenient derivatives to conduct our analysis.Now it is important to utilize this method to explore some nuances of points near the tricritical one.We will demonstrate how the method responds to the crossover effect.

B. Crossover phenomena
We begin by simulating the average eigenvalue as a function of T /T C .However, our focus now shifts to examining points near the tricritical point (TCP) to observe how the spectra of Wishart matrices behave when approaching this point alongside time series of magnetization simulated with (MC) simulations.By repeating our procedure, we initially investigate the issue in two dimensions to understand how the spectrum respond to the expected crossover phenomena in this model (refer to Figure 7).To accomplish this, we employed the values D = 1.9, 1.92, 1.9336, 1.9421, 1.9501, and 1.96582 (TCP).
We notice that the minimum becomes less pronounced and deformed as we approach the TCP.
However, it is interesting to note that even for points near the TCP, the method indicates the critical point albeit with reduced precision.Initially, the peak transforms into a shell, resembling a shoulder, and eventually, at the tricritical point, the minimum completely disappears.highlighting its divergence at T = T C .Interestingly, it is observed that the inflection point appears to be even more pronounced in three dimensions.
This indicates that the average, which localized the critical points well away from the TCP, strongly suffers the influence of this point, showing that the spectra of our correlation matrices precisely reflects what occurs with the thermodynamics of the model.Following this, we observe the dispersion of eigenvalues.We plot (∆λ) 2 as a function of T /T C for the same values of D previously used to study ⟨λ⟩.This result is presented in Figure 8.
In contrast to the behavior observed with ⟨λ⟩, the quantity (∆λ) 2 exhibits an inflection point very close to T = T C , even for points near the TCP; i.e., the variance senses the crossover but is not completely extinguished as with the simple average.However, at this precise juncture, a peak occurs at the TCP, which, upon closer examination, appears to shift as D approaches D T CP , i.e., we observe a migration of the maximum that will coincide at the critical temperature only exactly at TCP. Particularly intriguing is the observation that the migration of the maximum occurs with a decrease in its amplitude as D approaches D T CP .Exactly as occurred in the two-dimensional version of the model, the method indicates the critical point but loses precision as it approaches the TCP.However, we observe that for the 3D version, the minimum is more persistent, even at the TCP itself, since we do not observe a shoulder as obtained in the two-dimensional version.
Here, it is important to mention that in the three-dimensional version of the BC model, the short-time regime presents a logarithmic correction [20,21], which should suggest such different behavior.However, again, for points far from the TCP, the minimum of ⟨λ⟩ always occurs at .Average eigenvalue approaching the TCP in the two-dimensional BC model.We can observe that the shape of the curve is deformed as we approach the TCP on the critical line.method appears to reasonably respond even for points closer to the TCP.We can observe the inflection point up to just before the TCP, but we also notice a small deviation between the critical exact values and those determined by the method due to the crossover.At this precise TCP, there is a peak at the tricritical temperature that shifts from the previous points.Interestingly, at the TCP, we do not observe the inflection point in two dimensions.from the exact critical value.This further reinforces and suggests that we can use the inflection point of the spectral variance as a reliable indicator of critical phenomena (see Figure 10) far from TCP, but crossover effects can generate small deviations around TCP.
Crossover effects are observed in several works, and they play an important role in determining other quantities related to critical behavior, such as critical exponents.Here, we study their influence on the spectra of Wishart matrices built with time series of magnetization of the BC model.
Our results suggest that ⟨λ⟩ works very well for critical points outside the influence of the crossover, but it is not a good indicator of criticality near the TCP.In this case, we can make use of eigenvalue variance, which exhibits an inflection point at the critical temperature and responds reasonably well even when near the TCP, although it is also sensitive to crossover effects.
It is important to mention that MC simulations, whether in equilibrium or nonequilibrium, are generally sensitive to crossovers.For example, the dynamic exponent z, expected to be universal, is significantly influenced along the critical line in two dimensions [23,24], and even in the mean-field regime [22].
In Statistical Mechanics, the role of the maximum eigenvalue appears in many contexts, and an important question is whether they can also be used to quantify critical phenomena in the spectral method developed here.In other words, does the maximum eigenvalue of Wishart matrices respond to the critical behavior of the BC model?The answer is positive, and we will present the results in the next subsection.

C. Analyzing Extreme Statistics of Correlation Magnetization Matrices
The utilization of extreme values has been extensively investigated within formal contexts to characterize phase transitions in random matrices (see, for instance, [49,50]).Nevertheless, we posit that our approach holds promise for extension, leveraging similar principles computationally and efficiently to pinpoint critical points.Hence, this paper embarks on an exploration of extreme value statistics as indicators of critical points within the BC model, employing our correlation magnetization matrices.
Accordingly, for each matrix G * constructed, we extract its maximum eigenvalue and compute the average across multiple runs using the following formula:  and consider its behavior as a function of different temperatures for the BC model in both two and three dimensions.For such analysis, we choose D = 0, 0.5, 1.0, 1.75, 1.9, and 1.92 in two dimensions and D = 0, 1, 1.5, 2, 2.4, and 2.52513 for the three-dimensional version of the model.
The behavior of ⟨λ max ⟩as a function of T /T C is shown in Figures 11 and 12 respectively for the cases of the two and three-dimensional BC models.
We observe that in both situations, the critical point is identified by a notable inflection point.
Additionally, we show the first derivative in relation to T T C , simply described as , as a function of T T C as inset plots in these figures.Thus, we can also observe that the averaged maximum eigenvalue responds to the criticality of the system for the different critical points studied here for the BC model regardless of dimensionality.
This adds an additional parameter to our framework to identify criticality in spin systems that can be tested in other models.In the next section, we will conclude our analysis by showing that the method that uses short times (in this current contribution and in the previous ones [36,37] we used N M C = 300 steps) also works with small systems.Up to now, we have used L = 100 in two dimensions.We will demonstrate that this number can be further reduced.

D. Finite Size Scaling: Exploring Small Systems with Short Time Scales
The method's efficiency in saving computer time through the use of short time scales presents a particularly intriguing prospect.For instance, in this current study, we employed N M C = 300 steps.Therefore, to highlight the versatility of our method, we will explore another aspect: the system size.We have investigated this aspect in both two and three dimensions, demonstrating that systems can be studied effectively with even smaller sizes, yielding good estimates.
We deliberately selected only the case of D = 0 without loss of generality.The average eigenvalue is plotted as a function of T /T C for different sizes of the two-dimensional BC system (refer to Figure 13).Initially, we explore sizes ranging from L = 2 to L = 16, and subsequently extend to L = 20, 25, 30, 32, 64, 100, and 128.
We can observe an influence of the system size, where for small systems, a minimum at the exact T C is found.The inset plot illustrates that for L ≥ 32, the minimum at T = T C coincides.With L ≥ 64, there is excellent agreement.
For the three dimensional model, our investigation yields similar results.Encouragingly, we found consistent behavior, particularly noteworthy for L ≥ 16, where a distinct trend emerges: the average eigenvalue reaches a minimum precisely at T = T C .This observation implies the feasibility  of exploring intricate phenomena within compact systems.Thus, the potential for fruitful spectral analyses in modest-scale systems becomes increasingly evident.
As we wrap up this subsection, it's remarkable to note that in addition to employing shorttime evolution of magnetization (N M C = 300 steps), we can also leverage small systems to achieve robust results.Surprisingly, for L ≥ 64 in the two-dimensional BC model and L ≥ 16 in the threedimensional BC model, our method accurately identifies the critical temperature of the model.This underscores the efficacy of the spectral method, highlighting its strength in pinpointing critical points.
Readers are encouraged to juxtapose our method with simpler MC simulations, considering its additional workload in terms of matrix diagonalization.However, it operates on low-dimensional matrices (with N sample = 100 here -adjustable for further optimization).In comparison to timedependent simulations, which may require, for instance, L = 256 and a substantial number of runs to adequately sample quantities (with a minimum of 2000 runs for ferromagnetic initial states and over 10 4 for disordered initial states where m 0 ≈ 0), our method presents an intriguing alternative.
It's worth noting that capturing the thermodynamics of the model with such accuracy using a computationally "cheap" spectral method is not a trivial achievement.
Equilibrium MC simulations are plagued by the issue of critical slowing down, compounded by the use of larger lattices than those employed in our approach.While a thorough comparison between this spectral method and standard MC methods warrants attention, we must emphasize the compelling observation that the thermodynamics of the systems are remarkably well-reflected by this "spectral thermodynamics".
An avenue ripe for exploration is the investigation of long-range systems, which will undoubtedly command our focus in future applications, precisely due to the lack of requirement for large-scale systems.

IV. CONCLUSIONS
In this study, we have extended a method originally developed in [36] to describe the spin-1 Ising model with anisotropy, known as the Blume-Capel model.This model exhibits a tricritical point in both two and three dimensions.Our method has proven effective in accurately capturing these critical points and illustrating the associated crossover phenomena.
Furthermore, we underscore the computational efficiency of our proposed method compared to similar approaches.By diagonalizing matrices of size O(N sample ), where N sample is set to 100 in this work, we alleviate the computational burden.This stands in contrast to other spectral methods in the literature, which necessitate diagonalizing matrices of size O(L d ), where L represents the linear dimension of the system and d its dimensionality.For instance, in a system with L = 100 and d = 2, this would entail diagonalizing matrices of size 10 4 × 10 4 for a significant number of runs, which is computationally intensive.
In summary, our findings demonstrate that spectral methods provide a promising avenue for characterizing the thermodynamics of spin systems exhibiting tricritical points and crossover phenomena, regardless of the system's dimensionality.Notably, we achieved these results using very small systems and short time series, suggesting a means to bypass both critical slowing down and the necessity for extremely large systems often observed in standard MC simulations that do not involve the diagonalization of Wishart matrices.
While our proposal does not seek to directly compete with standard MC simulations, either in equilibrium or nonequilibrium settings, our results indicate that the method merits consideration for application in these contexts due to its efficiency and sensitivity.
Lastly, we emphasize the remarkable success of the developed model in characterizing chaos [51], as well as in describing the aging effects in spin systems [52], underscoring the ongoing exploration of its full potential and the depth of understanding yet to be achieved in this research domain.

V. APPENDIX
We observed that the inflection point appearing in the spectral variance serves as an excellent indicator of critical behavior when employing the method of constructing Wishart matrices.Clearly, within this framework, we can explore the nature of the inflection point under consideration.As a preliminary test, we opted for the simple case D = 0 and plotted the spectral variance for various values of t = T /T C , including t = 1, as depicted in Figure 15.We fitted two suggestive functions to the data, beginning with the well-known logistic function: and the second one, the Boltzmann function: Both fits were conducted by initially fixing t 0 = 1, as c 1 and c 2 represent normalization and fitting parameters, respectively.Thus, the fit essentially revolves around one parameter: p.This observation arises from the fact that both functions can be expressed as: where φ(t) = 1 1+t p for the logistic function and 1 1+exp t−1 p for the Boltzmann function.We can observe a good visual fit, both with a coefficient of determination r ≈ 0.998.For additional information regarding the fits presented in Figure 15, for the logistic function we obtained: c 1 = 50.4± 0.8, c 2 = 2.3 ± 0.4, and p = 11.8 ± 0.5.Similarly, for the Boltzmann function, we obtained: c 1 = 52 ± 1, c 2 = 2.8 ± 0.4, and p = 0.087 ± 0.004.
First, we would like to demonstrate that t 0 = 1 serves as an inflection point for both functions.
In the first case, we have: and it is certain that p > 1 for this type of curve.Thus, for t > 1, d 2 φ dt 2 > 0, while for t < 1, dt 2 < 0, demonstrating that the Boltzmann function with p > 1 has an inflection point at t = 1.Specifically, d 2 φ dt 2 t=1 = 0.
Similarly, for the logistic function, we have: = 0.
In the same Figure 15 we purposively presented an additional fit to a function: var(t, t 0 ) = c 1 (t − t 0 ) 1/3 + c 2 .
We tested other powers except for 1/3, as it provokes interesting discussions among more skeptical readers, even though it does not visually fit well.The aim here is to demonstrate that even when we encounter an inflection point where the second derivative is not defined, as suggested by the plots presented in Figure 6, we encountered no issues.For the best fit, we obtained c 1 ≈ −36 and c 2 ≈ 29 for t 0 = 1 (fixed).Despite d 2 var(t,1) dt 2 = 5 9 c 1 (t − 1) −5/3 and t 0 = 1, ensuring an inflection point even when:

Figure 1 .
Figure1.The phase diagrams for the two and three-dimensional BC models are depicted.The points utilized in our numerical experiments are extracted from Butera and Pernici[43], serving as foundational data for the investigations conducted in this study.
Following this, in the second subsection, we extend our investigation to demonstrate the crossover effects in the model, as captured by the density of maximal eigenvalues of Wishart matrices.For our analysis, we construct N run = 1000 distinct matrices G * of size N sample × N sample for each fixed temperature.Each matrix is derived from N sample = 100 magnetization time series, each comprising N M C = 300 Monte Carlo steps.These time series are obtained via MC simulations employing heatbath single spin flip dynamics for the BC model, resulting in a total of 10 5 eigenvalues used to construct the histogram for each temperature.All eigenvalues are categorized into N b = 100 bins.In the two-dimensional systems, we utilize a linear dimension of L = 100, while in three dimensions, we employ L = 22.

Figure 2 .
Figure 2. The Density of States in the Two-Dimensional BC Model with anisotropy (D = 1).The gap between eigenvalues varies with the temperature of the simulated system.While the system approaches the MP law, an exact match is not achieved at high temperatures (T > T C ) due to the presence of spin-spin correlations, preventing complete correspondence.

Figure 3 .
Figure 3.The density of states for anisotropy D = 1 in the three-dimensional BC model.Similar behavior

Figure 4 .
Figure 4. Average and variance of the two-dimensional BC model as a function of temperature are depicted.The inset plots show the derivative of the variance, indicating a divergence at T = T C .

Figure 5 .
Figure 5.The average and variance of the three-dimensional BC model as a function of temperature illustrate a similar behavior occurring in three dimensions.The inset plots depict the derivative of the variance,

Figure 6 .
Figure 6.Second derivative of variance (ζ) for both the two-dimensional and three-dimensional BC models.The critical temperature precisely corresponds to the inflection point of the eigenvalue variance, indicated by ζ < 0 for T < T C and ζ > 0 for T > T C .

Figure 7
Figure 7. Average eigenvalue approaching the TCP in the two-dimensional BC model.We can observe that

Figure 8 .
Figure 8. Eigenvalue variance as a function of temperature approaches the TCP in the 2D-BC model.The

Figure 9 .
Figure 9. Average eigenvalue approaching the TCP in the 3D BC model, mirroring the analysis conducted for the 2D version shown in Figure 7. max

Figure 10 .
Figure 10.Eigenvalue variance approaching the TCP in the 3D-BC model.We can observe the inflection point until the TCP point, but that slightly differs from the best estimates of the critical temperatures in this vicinity of the TCP.

Figure 11 .
Figure 11.Averaged maximum eigenvalue as a function of T /T C for different values of D in the twodimensional BC model.

Figure 12 .Figure 13 .
Figure 12.Averaged maximum eigenvalue as function of T /T C for different values of D in the threedimensional BC model.

4 ⟨𝜆⟩Figure 14 .
Figure 14.Average eigenvalue plotted against T /T C for various linear system sizes.We consider L = 2, 4, 8, 10, 16, and 22 in a three-dimensional BC model with D = 0 for simplicity.The inset plot highlights that for L ≥ 16, the minimum occurs precisely at T = T C .

Figure 15 .
Figure 15.Spectral variance as a function of various values of t = T /T C for D = 0, including t = 1, is presented.We provide fits using the Boltzmann, logistic, and a function exhibiting a unique type of inflection point at t = 1.
< 1, and d 2 φ dt 2 < 0. This shows that t = 1 is also an inflection point.Again d 2 φ This scenario is a basic calculus concept, yet it is always worth noting.